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Abstract: We calculate the deconfinement line of transitions for large Nc QCD at 
finite temperature and chemical potential in two different regimes: weak coupling in 



O ■ the continuum, and, strong coupling on the lattice, working in the limit where Nf is of 

■ order Nc- In the first regime we extend previous weak-coupling results from one-loop 

■ perturbation theory on S*^ x S*^, where the theory reduces to a matrix model, analogous 



to that of Gross, Witten, and Wadia, to higher temperatures. We obtain the line of 
transitions that extends from the temperature-axis, where, to a first approximation 
the transition is higher than fourth order, to the chemical potential-axis, where the 
^ ■ transition is third order. In the second regime we use the same matrix model to 

■ obtain the deconfinement line of transitions as a function of the coupling strength 

and fi/T to leading order in a strong coupling expansion of lattice QCD with heavy 
quarks, extending previous U{N^ results to SU{N^. We show that in the case of zero 
chemical potential the result obtained for the Polyakov line from QCD on S"^ x S'^ 
at weak coupling, reproduces the known results from the strong coupling expansion, 
under a simple change of parameters, which is valid for sufficiently low temperatures 
and chemical potentials. 



1 Introduction 



Obtaining the phase diagram of QCD at non-zero chemical potential is a long-standing 
problem, and, although complicated, results are starting to emerge which suggest that 
there is a range of new phenomena at non-zero density that demand new explana- 
tions. The presence of the sign problem, however, in combination with a large coupling 
strength at low and moderate temperatures and chemical potentials, has the conse- 
quence that conventional techniques such as lattice simulations based on importance 
sampling, and infinite volume perturbation theory, are only able to access a small frac- 
tion of the overall phase diagram. To obtain a better understanding of the difficulties 
faced at non-zero chemical potential, and to obtain a picture of the phase diagram 
where conventional techniques break down, it is helpful to study limits of QCD which 
make calculations analytically tractable. One such limit is obtained by taking A^^ large, 
which introduces several benefits: large factorization simplifies the action by remov- 
ing terms that involve correlations between different points in spacetime, the infinite 
number of degrees of freedom provided at large Nc allows for sharp phase transitions 
in finite volumes, and the deconfinement phase transition, which takes places at suffi- 
ciently high temperature, becomes describable in terms of the behavior of a continuous 
distribution of eigenvalues of the Polyakov line order parameter, which reduces the 
theory to an analytically solvable matrix model. 

In what follows we describe in detail a calculation to obtain the Polyakov line 
observable and map out the confined and deconfined regions of the large N^. QCD 
phase diagram. First, we consider QCD on x S*^, where the radius, i?, of the 5''^ is 
taken sufficiently small that perturbation theory is valid, R <C ^qcd- ■'■^ comparison 
with S*^ X M^, the infinite spatial volume has been sacrificed in exchange for being able 
to work at any temperature, T, where T = ^ and (3 is the length of the S^. In [1] the 
phase diagram for QCD with chemical potential on x was calculated from various 
observables in the low temperature limit. In the first part of this work, we extend these 
results to higher temperatures to obtain the deconfinement transition line in the fi-T 
plane from the /i = axis to the T = axis, where we match onto the low temperature 
result in [1]. 

The differences in the weak coupling analysis in this work from [1] are 1) we keep 
the finite temperature contributions in the action, but 2) we will make the approxima- 
tion that z^n = 0, and 2j„ = 0, for n > 1, where z^n and Zfn are the single particle 
partition functions for vectors and fermions, for a given number of windings n. This 
approximation was also used in sections of [2] to obtain concrete results at finite tem- 
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perature where the corrections are exponentially suppressed at low temperatures, and 
power law suppressed at high temperatures. In taking this approximation we restrict 
our weak-coupling results to the region of temperatures which are not too high, and we 
restrict the chemical potential to /i ^ ^/iLi?~o ~ ^^^^^ ^fi ^ + + m?R^ 

are the energy levels of fermions of mass m formulated on S*^, where we consider the 
case where mR ~ 0. This includes the region of the deconfinement line of transitions 
extending from the T axis to the yU axis for the case of very light quarks. To obtain re- 
sults for larger chemical potentials it would be necessary to include terms with z j„ 7^ 
for larger n. To consider larger quark masses it would be necessary to include terms 
with 7^ for larger n. We investigate how such corrections would be included in 
a manner similar to that of [2] and give some results for Yang-Mills theory in Section 
3.2.1. 

To obtain finite temperature results we follow [1] and generalize the technique of 
Gross and Witten [3], and Wadia [4], to include the contribution from the quarks cou- 
pled to non-zero chemical potential. We start by reviewing Yang-Mills theory in Section 
3, and QCD with zero chemical potential in section 4. In each case we derive the weak 
coupling result on S"^ x S'^ considering a more general scenario where the eigenvalues 
of the Polyakov line are allowed to lie off of the unit circle until the end. Then we 
use a simple change of parameters to show that the matrix model also reproduces the 
Polyakov lines for Yang-Mills theory and QCD with /i = from the strong coupling 
expansion with heavy quarks in [5]. The calculations for /i 7^ are given in Section 5. 
An important feature, also imposed in [1], is that the gauge fields are complexified in 
the case of 7^ 0. Thus we allow the Polyakov line eigenvalues to lie off of the unit cir- 
cle throughout. This is crucial in obtaining the correct stationary point solution since 
the action of QCD with /z 7^ is complex (This is also important in lattice simulations 
at finite density using stochastic quantization and complex Langevin methods. See, for 
example [6, 7]). It is also critical to impose the SU{N^ constraint to obtain a solution 
which results in a /^-dependent free energy and a non-zero quark number, as done in 
[1] for QCD on 5*^ x S*^, and in [8] for effective spin models treated with mean field 
theory. For the U{Nc) theory, at least for QCD on x S^, in effective spin models, 
and in the lattice strong coupling approximation [9], the non-zero chemical potential 
results in a trivial shift of the gauge field, resulting in Polyakov lines which are given 
by p±i(/i)=p±i(/i = 0)e^'^'5. 

After obtaining the Polyakov lines at weak coupling on 5*^ x S*^ as a function of 
the temperature and chemical potential and calculating the phase diagram we use the 
same matrix model to obtain the Polyakov lines from a lattice strong coupling expansion 
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with heavy quarks as a function of the couphng strength and the ratio yu/T, in Section 
6. Our resuhs extend the results in [5, 9] to SU{Nc) QCD with 7^ by imposing 
the corresponding Ylf=i = constraint on the eigenvalue angles of the Polyakov line. 
This results in a non-zero quark number, and a /^-dependent free energy. The change of 
parameters which makes conversion between weak and strong coupling results possible 
at sufficiently low temperatures and chemical potentials is derived next in Section 2. 
It is valid for observables such as the Polyakov lines, quark number, and free energy, 
which don't include correlations between different spacetime locations, in other cases 
the procedure would need to be generalized, if possible. It is unclear to us how it comes 
about that such an approximation exists which leads to the equivalence of Polyakov 
line matrix models connecting observables in weakly-coupled, small volume continuum 
QCD, to strongly-coupled large volume lattice QCD. It is also not clear, and would be 
interesting to find out, if such an approximation is also possible for other theories. 



2 1-loop QCD on x vs. lattice strong coupling expansion 
with heavy quarks, as Nc ^ 00 

The action of large Nc QCD, with ^ fixed, to leading order in the lattice strong 
coupling expansion and the hopping (heavy quark) expansion is given by [5, 9, 10] 

Slat - Svdm = - JDJ2[{W)W\x) + {W^)W{x) - {W){W^)] 

(2.1) 

X 

(a \Nr ^ 

where Svdm is the Vandermonde contribution to the action, J = 2 ( j in terms of 

N 

spacing, and D is the number of spatial dimensions, and W{x) = Tr J^^Jq"^ Ut^i is the 
Polyakov line. 



the inverse coupling Piat = and the number of lattice sites in the temporal direction 
Nt, h = 2^K^^ is the hopping parameter with k = , where a is the lattice 



On X the action from one loop perturbation theory takes the form [1, 2] 

00 ^ 

Ss'^xS'-^ ~ Svdm = — N —Zynpnp-n 



n 

00 (2-2) 
+ ^^^^ E H^^^" + , 
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where 



1=1 



'fn 



(2.3) 



2^/(/ + l)e-"4V('+i)^+™^^% (2.4) 



and the Polyakov hnes are defined by p„ = ±-Ti ^ = ^^e'"^" = j^- e^"^% 
where a = ^ J^^g d^xAo(x) and dta(t) = 0, such that a gauge with constant, diagonal Aq 
is chosen and massive, off-diagonal fluctuations have been integrated out. In (2.2) the 
Yang-Mills contribution has pnp-n in place of {W)W^{x) + W{x){W^) - {W){W''). In 
practice the equations of motion obtained from derivatives on {W)W'' {x) + W {x) {W'^) — 
{W){W'^), with the expectation values flxed, are equivalent to those obtained by taking 
the derivatives on pnP-n- The actions (2.1) and (2.2) are otherwise similar with the 
exception of the sum over n in the formulation on x S^, and that the x-dependence 
is kept in the strong coupling expansion. Due to the absence of terms with correlations 
between different sites in the action (2.1) an observable of the form {F(W, W"^)) can be 
obtained as follows 

{F{W,W^)) = J-^ I JJdiy(x)e-^[^(")'^'(")]^F[iy(x'),iv^(x')], 

_ /diye-^(^'^^)F(iy,iyt) ^ • ^ 

~ / dVTe-'^Ww/t) 

Therefore, when it is possible to approximate the sum over n in (2.2) by the n = 1 
contribution, then it is possible to calculate observables of the form (F(pi,p_i)) in 
weakly-coupled QCD on x 5''^ and use the transformations 

Nc^ " (2.6) 
z^i ^ JD , 



'/I 



^ h. 



to obtain the result for strongly-coupled lattice QCD with heavy quarks, or vice versa. 
The approximation of keeping only the n = 1 contribution corresponds to taking 
Zvn = 0, Zfn = 0, for u > 1, which is a good approximation for p < eji and when 



- 5 - 



the temperature is not too high ^. Whether the corrections, corresponding to terms in- 
cluding multiple windings of the Polyakov lines, can be included by considering higher 
order terms in the strong coupling and hopping parameter expansions such that the 
weakly-coupled and strongly-coupled theories still share the same matrix model is left 
for future research. 



3 Yang-Mills theory 

Taking Nf = and evaluating the Vandermonde contribution in (2.2) leads to the 
one-loop action for Yang-Mills theory on x [2] 

OO 

N^y^ -{1- Z,n)pnp-n. (3.1) 



Sym 

n=l 



Taking the Nc ^ oo limit makes it possible to define a map 



1 

— T 



where the Polyakov line eigenvalues Zj = e are distributed along the contour C with 
a density 

q{z)^z^. (3.3) 
d^ 

From the map in (3.2) it is clear that the density g{z) is complex but evaluated along 
the contour C the quantity ^dzg{z) is real and positive and normalized as 



(3.4) 



so, as n increases from the Nc eigenvalues along the circle move off into the complex 
plane but still lie on a 1-dimensional arc, as we expect from the general behavior of 
matrix models with complex potentials, for example [12]. The map can be used to 
simplify the equation of motion obtained from ^ = to the form 



^ / S^'(^')7^ = E^-(/^-"^"-/^"^"") ' (3-5) 

n=l 

where ^ indicates that the principal value is taken with the point z left out of the 
range of integration, and the contour C is an arc which opens on the negative real-axis 
with endpoints at z = re^'^ , z* = re~^^. 



"'^Notice that in this limit the matrix model reduces to the one in [11], for which the effective 
potential and Polyakov lines were obtained for Nc — 3. 
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3.1 Confined (ungapped) phase 

Following [3, 4] we represent the confined phase with a continuous, ungapped dis- 
tribution of the Polyakov line eigenvalues, such that the contour C, along which the 
eigenvalues are distributed, is closed and the endpoints z, z* coincide with ip = n. The 
equation of motion (3.5) can be evaluated using Cauchy's theorem. Fourier expanding 
the density as 



and collecting the residues at z' = 0, z, results in p„ = p_„ = for n 7^ 0. This is 
consistent with the result {W) = within the confined phase from the strong coupling 
expansion of Yang- Mills theory. Using the identity (3.4), it is clear that po = 1. From 
the weak coupling action (3.1) the confined phase persists while z^i < 1 [2], or, using 
the transformations in (2.6), the confined phase persists while JD < 1 in the case of 
strong coupling, as found in [5]. 

3.2 Deconfined (gapped) phase 

Following [3, 4] the deconfined phase is obtained when the distribution of the Polyakov 
line eigenvalues develops a gap on the negative real axis, such that C is open with 
costp 7^—1. We review the result in [2] for the Polyakov line from Yang-Mills theory 
on X S^, in the deconfined phase, leaving open the possibility that the distribution 
of the Polyakov line eigenvalues lies off the unit circle, and that pi 7^ p_i, until the end. 
The singular integral on the LHS of the equation of motion (3.5) is solved by defining 
a resolvent 



00 




(3.6) 



n=l 




(3.7) 



Using the identity (3.4), the resolvent can be rewritten as 




(3.8) 



Therefore 0(2;) satisfies the Plemelj formulae 



(P+iz)-(P~{z) = 2zg{z), 



(3.9) 



00 




(3.10) 



by taking the contour C to lie on a square root branch cut. Following [4] the Plemelj 
formulae are solved writing 



with 



From (3.10) 



[z) = h{z)H{z), (3.11) 



h{z) = ^/ {z - z){z - z*) . (3.12) 



2mf^z'-z ^^z' -~z){z' -~z*) ^ ' 

where F is a contour around C which can be peeled off to evaluate the integral by 
collecting the residues at 0, oo, and z to obtain 



H{z) = [{z -z){z- z*)]-''^ J2 (Z^--^" - P-^"") 

(3.14) 



1=1 k=0 

From (3.11) this results in 



oo 
n=l 



^ Zyn {p-nZ^ - pnZ ") + \/ {z - z){z - Z*) 

oo oo (3.15) 

X ^^Pfc(cos?/^)2„(i+fc) (p/+fcr"''"^z"' + p_;-fcrV-^) , 



1=1 k=0 

which agrees with [13], where Pk{x) are the Legendre Polynomials. The density is 
obtained from (3.9) as 

g{z) =a/ {z - z){z - z*) 

X ^ ^ Pfc(cos V)2^a+fc) {pi+kr-''~^z-^-^ + p_i_kr^z^-^) . 

1=1 fc=0 



The Polyakov lines can be calculated from 



dz 

Pn= / ^f^(^)^", (3.17) 



c 
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where again we write the integration along the C as the contour F, then peel T off and 
evaluate the integral collecting the residues outside. Integrating both sides of (3.9) one 
obtains the more general transformation 

such that the Polyakov lines can be obtained from 

Pn = f^^/{z)z\ (3.19) 

where po is the identity (3.4). If we make the approximation that z^n = for n > 1, 
then 

Po = 1 = ^z^i (pir"^ + p_ir) (1 - cosip) , (3.20) 

gives the constraint 

pir-i + p_ir = -. (3.21) 

^1,1(1 — costp) 

Similarly one can solve (3.19) for = ±1 to obtain 

ir^2,,ip_i sin^-^ 
= ^Hli 1- , (3.22) 



2 — z^i{l — costjj) 
\r'^ZviPi sin^?/' 

2 — Zyl{\ — COS'?/') 



P-i = n M TT ' (3.23) 



which imply the constraint 



2 — z„i(l — cos '?/') = -2^,1 sin^ . (3.24) 
The solution from this constraint is 

cosV' = -1 + — Jzl-^ - z^i , (3.25) 

such that (3.21) gives an expression for the Polyakov lines as a function of the temper- 
ature 



1 



Pir-^ + P-ir = 1 + Jzl^-z^^. (3.26) 
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Figure 1. Polyakov line pi as a function of temperature for Nc = oo Yang-Mills theory on 
X S"^. The red line corresponds to the result from the confined phase. The blue curve 
corresponds to the result from the deconfined phase. 

For Yang-Mills theory the action is real so we take r = 1 and p_i = pi such that the 
Polyakov lines simplify to 



which is the result in [2] and we plot it in Figure 1 along with the p±i = result 
for the confined phase from the previous subsection to illustrate the discontinuity at 
the transition temperature. The transition occurs when z^i = 1 which corresponds to 
TR ^ 0.759 where p±i jumps from to |. As TR — )• oo, z^i — )■ oo and p±i — )■ 1. 

Notice that under the change of parameters (2.6), (3.27) becomes the known result 
for Yang-Mills theory from the lattice strong coupling expansion [5] 



where the transition point z^i = 1 is converted to JD = 1, as in [5]. 
3.2.1 Corrections 

It is straightforward to consider the Zy2 7^ corrections to the Polyakov lines for Yang- 
Mills theory in the deconfined phase. In this section we simplify our notation by taking 







0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 
TR 




(3.27) 




(3.28) 
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X = costp. The identity constraint given by (3.19) for n = becomes 
1 = ^(1 - x) [2z^i{pir~^ + p_ir) + 2:^2(1 + 3x)(p2r"^ + P-2r' 



(3.29) 



Because the action is real we take r = 1 which leads to p„ = p_„,. (3.19) for n = 1, 2 
gives 

4(1 + x)2 



Pi 



P2 



4(l + 3a;) + z,i(l-x)3 
8-2z^i(l -x)(3 + x) 



(3.30) 
(3.31) 

(3.32) 



(l-x)z,2[4(l + 3a;) + z,i(l-x)3] ' 
where x = cosip is the solution of 

^ ^ 8(l-x)^(l + a:)"z,iz,2 

[4 - z^i(l - x)(3 + x)] [16 + 2^2(1 - x) [7 + X (3 + x(13 + 9x))]] 

which is real and has |x| < 1. The correct root depends on ^: it changes when z„2 = 1, 
corresponding to TR ^ 1.519. However, the Polyakov lines pi, p2 remain continuous 
over TR ^ 1.519, as shown in Figure 2. Notice that pi in Figure 2, which contains the 
Zv2 7^ corrections, has only barely changed from the result with z^2 = in Figure 1. 



4 QCD with 11 = 

The one-loop action for QCD with p = is [2] 

SqCD = ^ - (1 - Z^n) PnP^n + NfN^ ^ Z/„ (p„ + p-„) 



n=l 



n=l 



such that the equation of motion becomes 
dz' I z' + z 



I ^vnPn 



^'-lTzf„,\z-^ 



(4.1) 



(4.2) 
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Figure 2. Polyakov lines pi, p2 as a function of temperature for Nc = oo Yang-Mills theory 
(including contributions resulting from z^n 7^ for n = 1, 2.). The blue curve corresponds 
to pi (note that the there is also a disconnected section of the blue curve along the TR axis 
until z^i = 1). The red curve corresponds to p2- 

4.1 Confined phase 



As in Yang-Mills theory the contour C is closed in the confined phase and the equation 
of motion (4.2) is evaluated by Fourier expanding the density as in (3.6) and collecting 
the residues at 0, z, resulting in 



P±n 



■fn 



Nc 1- z, 



(4.3) 



for n 7^ 0, in agreement with [14]. For n = 0, po = 1 as in Yang-Mills theory. For n = 1 
we can obtain the lattice strong coupling result using the transformations in (2.6) to 
obtain 



■{W) 



h 



1 - JD 



(4.4) 



in agreement with [5]. 



There is an upper limit on the Polyakov line that can be obtained while the dis- 
tribution remains ungapped. A consequence of working in the approximation where 
= Zfn = for n > 1, is that the Polyakov lines, Pn = for n > 1 in the confined 
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phase. The eigenvalue density (3.6) obtains a zero, 



q{z) = p^i + + piz-^ = Q , (4.5) 

when 

-1 ± VI - 4pip_i 



2p-i 



(4.6) 



hes on the Polyakov hne eigenvalue distribution, on C. Since the distribution lies on the 
unit circle when /i = 0, with p_„ = p„, the gap forms at z = —1, which corresponds to 
Pi = p_i = |. This is the latest point at which the ungapped phase can end. Plugging 
P-ti = ^ into (4.3) gives the allowed region of the ungapped phase as 

l-z,,-2zj,j^>0. (4.7) 

4.2 Deconfined phase 

After a gap forms in the Polyakov line eigenvalue distribution the theory enters the 
deconfined phase, where the contour C lies on an arc which opens up on the negative 
real-axis. Comparison of the equation of motion (4.2) with the equation of motion 
for Yang-Mills theory (3.5) indicates that there are no new poles in z so the results 
obtained in the deconfined phase for Yang-Mills theory should carry over to QCD with 
p = by shifting z^np±n z^np±n - -j^i-lTzfn. 

Taking z„„ = z/„ = for n > 1, and defining x = cosip, the constraint from the 
identity, (3.19) with n = 0, becomes 



Po = l= ^(1 



Nf. , f. . , Nf. 



(4.8) 



For n = ±1 (3.19) gives the Polyakov lines 



_ ;z/i^(l - x) [8 + 4r±2(i + x) - z,,{l - x)2(3 + x)] 
~ 16 - z,i(l - x) [16 - z,i(l - x)2(3 + x)] 



(4.9) 



Since p = the action is real and the eigenvalues of the Polyakov lines lie on the unit 
circle with r = 1. In this case 



Pi = P-i 



(1 -x)(3 + x)z/i§ 



4- {l-x){3 + x)z^i 



(4.10) 
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Solving the constraint (4.8) for x gives 



_ -^vi-2zfi^ + 2^Jzl^- Zyi + 2zyiZfi^ + zj^^ (4.11) 

X , 



where we have chosen the root that gives |x| < 1 as a function of T. Plugging this into 
(4.10), allows us to solve for p±i, 



Zvi zfi -t- Y ^1,1 ^vi -T ■^^vi^fi -i- ^/i (4-12) 



ZZyl 



which is the result in [14]. This reduces to the Yang- Mills result (3.27) when ^ — )■ 0. 
It also matches onto the result in [5] for yU = under the transformations (2.6), that is 



The equation for x{T) in (4.11) only has a solution with |x| < 1 for 

1-2,1-22/1^ <0. (4.14) 

The minimum temperature occurs when the LHS=0. From (4.7) this is precisely the 
temperature at which the confined phase must end and plugging it into (4.12) results 
in p±i = 1. Therefore the transition is smooth and occurs at precisely 

l-2,„i- 22/1^ = 0, (4.15) 

as in [14, 15], or in the case of the strong coupling expansion, using the transformations 
in (2.6), the transition occurs when 

l-JL)-2/i = 0, (4.16) 

in agreement with [5]. 

The weak-coupling result for the Polyakov line (4.12) is plotted in Figure 3 for 
mR = along with the result in confined phase from the previous subsection (4.3) 
where the transition point occurs at the temperature given by (4.15), corresponding to 
Pi = 1. This transition is at least fifth order, and could potentially be infinite order, 
resulting in a crossover, since a pattern appears to emerge in the derivatives of the 
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Pi 




0.6 0.8 
TR 



Figure 3. Polyakov line pi as a function of temperature for QCD with Nc, Nf = oo, = 1, 
and mR = 0. The red curve corresponds to the ungapped phase. The blue curve corresponds 
to the gapped phase. 



partition function with respect to the temperature, evaluated at the transition point ^. 
We have found that 



1 d''{\ogZ) 



2N^ ' 



(4.17) 



where k refers to the number of derivatives with respect to T, and we have checked that 
this is true up to /c = 4. Here, logZ is obtained from (4.1) in the hmit z^n = Zfn = 
for n > 1, and using S = — logZ, valid at the saddle point. (4.17) holds in both the 
confined and deconfined phases, where logZ is evaluated with pi obtained from (4.3) 
and (4.12), respectively. The resulting ^ g^k^^ in each phase is evaluated transition 
point (4.15), satisfying the derivatives of (4.15) as well. For the lattice theory (4.17) 
becomes 



1 9^(logZ) _ 1 
iVf dT>' ~ 2 



(4.18) 



It is unclear how the order of the transition would be affected when corrections are 
taken into account from z^n and 2/„ for n > 1 for the theory on x S^, or additional 
terms in the strong coupling and hopping parameter expansions for the lattice theory. 
This is a possible direction for future research. 



^We would like to thank Andras Patkos for bringing this possibility to our attention. 
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It would be interesting to see what happens to the order of the transition as the 
quark mass is varied. As m — )■ oo the quark contribution becomes less significant 
and the theory approaches pure Yang-Mills. To accurately determine the order of the 
transition for the theory on x for larger quark masses, whether the transition 
remains continuous or not to high order, it would be necessary to include corrections 
from z^n 7^ for n > 1, which, one by one, become of the same order of Zfi when the 
mass is increased. In the strongly coupled lattice theory this corresponds to including 
corrections to the leading order in the hopping expansion. 

We know from the low temperature results for /i ^ e/i in [1] that the transition 
becomes third order in the limit fi Efi. In what follows we expect a sharpening of 
the transition for the theory on x S''^ as a function of TR for increasing /iR, or in the 
lattice theory as a function of J for increasing /i/T, but it is difficult to obtain the order 
of the transition between the /i = 0, and T or J = endpoints, because the derivatives 
of log Z become complicated and need to be solved numerically. We are only be able 
to determine that the order of the line of transitions in the region < < Syi is at 
least third order. 



5 QCD with Id < £fi 

The one-loop action for QCD with /i 7^ is [1] 

00 ^ 

SqCD =NI - (1 - Z^n) PnP~n 



n=l ^ 

n=l i=l 

where AT is a Lagrange multiplier which imposes the SU{N^) constraint ^'.^^ Oi = 0. 
The equation of motion becomes 



•^'^ n=l 



(5.2) 



where we have defined a±n = ZynP±n — jfi — ^)^Zfne^"^^. 
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5.1 Confined phase 



In the region of the confined phase with n < Sfi solving the equation of motion (5.2) 
with the density (3.6) gives 

P±n = if ^ \ ^" . (5.3) 

Nc 1 - Z^n 

in agreement with [1]. Considering n = 1 and making the replacements in (2.6) gives 
the results for the lattice strong coupling theory 

J_(W) = — , (5.4) 

^(W^\ = J}^ (5.5) 

in agreement with [9]. 
5.2 Deconfined phase 

The lack of new poles in the equation of motion (5.2) resulting from the contribution 
of quarks with fi ^ means that the resolvent takes the form in (3.15) with the shift 
ZvnP±n — ^ a±n, such that it becomes 



n=l 

oo oo 



oo oo 

X 

1=1 k=0 



Using (3.19) with n = and working in the limit = ;z/n = for n > 1 we obtain 

po = 1 = ^ (1 - x) (air^^ + a_ir) . (5.7) 



The Polyakov lines obtained from (3.19) with n = ±1 are 

Zf 1^(1 - x)e^t'l^ [8 + 4e±2M/3^±2(^ + x) - 2„i(l - xf{'i + xjj ^ 



16 - ;2^i(l - a;) [16 - z^i(l - s)2(3 + a;)] 
Plugging (5.8) into (5.7) results in a relationship between r and x given by 

2z^i^(l-a;)(l + r2e2'^^) 
~ [4 - z^i{l - a;)(3 + x)] ' 



(5.9) 
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The solution for r which allows for matching onto the /i = result is given by the 
positive root, 



^ _ 4 - z,^{l - x)(3 + x) + ^[4 - z,i(l - x)(3 + x)f - 16(1 - x^z),^^ 

Another two constraints are necessary to obtain x or r as a function of the temperature 
(two because there is also a dependence of x and r on M). One constraint is that the 
solution must satisfy [13, 16] 

( Er=l(«-n^"-«n^'")-A/"+l 

which can be evaluated to obtain 

oo oo 

^ afc+iP,(x)r-('=+^) = 1-Af+J2 o^-kPk{x)r^ • (5.12) 

A:=0 k=l 

This agrees with the form in [2, 13] for A/" = 0, Nf = and r = 1. At least for the 
concerned case with = Zfn = for n > 1 this constraint is equivalent to that which 
results from expanding the resolvent in the z ^ oo limit. The 2; — )■ and z — > cxd 
expansions can be performed by using (3.7), (3.4), and (3.17), and these take the form 

00 

lim0(2) = 1 + 2^2>_„, (5.13) 

n=l 

00 



lim (h(z) 



z^oo ' Z 

n=l 



The terms in (5.13), (5.14) can be matched against the corresponding terms of the 
2; — 7- 0, 00 expansions of (5.6) to obtain the constraints 

1 = —J\f — xair^-^ + a^ir , (5.15) 

— l = —Af — air^'^ + xa-ir. (5.16) 

Notice that (5.16) gives the same constraint as (5.12). Also the z^^ constraints from 
(5.13), (5.14) reproduce the Polyakov lines in (5.8). The constraint (5.15) could also be 
obtained using the identity constraint (5.7) with (5.16). Using (5.8) in (5.15), (5.16) 
and solving for N' gives 



^^ = (l + x) 



1 ^ 4z/i^r-ie-^^ [z,i - z.^x^ + 2r^e^''^{2 - z,i(l - x))] 



X 16 — z^i(l — x) [16 — z^i(l — x)^(3 + x)] 

(5.17) 



- 18 - 



where the equation for x{TR) obtained by plugging r from (5.10) into (5.17) reduces to 
that of QCD with = 0, if A/" = 0, as expected. For /i 7^ it becomes possible to have 
Af because it corresponds to the normalized quark number J\f = -^Nq as in [1]. 
This can be seen by summing (5.2) over all the eigenvalues, Zj = e*^^ , of the Polyakov 
line, and solving for the Lagrange multiplier to obtain 

^ Nc 00 



1=1 n=l 

Nc 00 



^J](-l)"^,„ (e 



c i=i n=l 



J. J Nc 00 



(5.18) 



1=1 1 = 1 

1 laiogZ 



1 



Since (5.11) - (5.16) have only introduced one new constraint it is necessary to impose 
an additional constraint to solve for A/", the SU{Nc) constraint. In the large Nc limit 
Ylf=i = takes the form 



^^(z) log(2;) = 0. 



(5.19) 



Defining the contour F around C and peeling it off to enclose the poles at and 00, 
and the branch cut of log z on the negative real axis, gives 



dz 

( 

2'Kiz 



[z) log(2;) = lim 



€->0 

r] — ^00 



dz 
2'Kiz 



[z) \og{z) 



dz 
( 

2'Kiz 



[z) \og{z) + 



d^ 



(5.20) 

which is similar to [1] except that there is no pole resulting from taking the sum over 
n in the equation of motion since we have truncated it. Plugging in (f){z) from (5.6) 
results in an expression which constrains the Polyakov lines and other observables to 
satisfy the SU{Nc) constraint, 

(1 + x) logr 



[l-x)-{l+x) log(^) 



(5.21) 



which is clearly for r = 1. It is now possible to solve for A/", r, and x as a function 
of TR. Plugging r from (5.10) into (5.17) and (5.21) gives two equations for as a 
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Figure 4. Polyakov lines pi and p-i as a function of the temperature for mR = 0, jj- = 1. 
(Left) p = 0.6/R, (Riglit) p = 1.2/R. 



function of x, which we refer to Mi, and M, respectively. Setting A/i = M gives the 
equation for x(T), which can be used to solve for the Polyakov lines (5.8). 

In Figure 4, the Polyakov lines p±i are plotted as a function of the temperature 
for fi = 0.6/R, fi = 1.2/R. The red part of the curve corresponds to the ungapped 
distribution, and the blue to the gapped distribution. The transition is indicated by 
the dotted line. The sharper transition of pi compared to p_i was also observed in 
the low temperature results in [1], where the fact that pi ^ p*^i for 7^ was also 
observed and the extent of the difference served as an indication of the severity of the 
sign problem. 

As the temperature is increased there is a certain critical temperature, Ti^^ at which 
it becomes possible for the distribution of Polyakov line eigenvalues to develop a gap. 
This occurs at the temperature at which x{T) develops a solution, in the range [—1, 1], 
which is not —1. However, up to a second critical temperature, Tc , the ungapped 
solution is also still possible. This second critical temperature is given by the solution 
of (4.15), which is /i-independent, which happens because the transition occurs from 
a phase in which A/" = where the Polyakov line can be shifted to absorb the /i- 
dependence. Therefore, to determine the temperature at which the gapped solution 
takes over it is necessary to compare the free energies for the gapped and ungapped 
distributions for the range of temperatures, Tc^^ < T < Tc^\ in which both are possible. 
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Figure 5. Free energy of the gapped and continuous distributions for the range of temper- 
atures Tc^^ < T < Tc^^ in which both are possible, for mR = 0, ^ = 1. (Left) // = 0.6/i?. 
(Right) /i = 1.2/R. 



This is done in Figure 5 for /i = 0.6/-R and ^ = 1.2/R, using 

(5.22) 



7 = —F = —TS ~ T 



Mr 



It is clear from Figure 5 that for both /i = 0.6/i? and /i = 1.2/R the gapped solution 
results in a lower free energy for the entire range of temperatures in which both the 
gapped and ungapped solutions are possible. Therefore as the temperature is increased 
the gapped solution will take over as soon as it exists. We have compared the free 
energies for several other values of /i and this seems to hold for the full range = to 

/i ^ Sfi. 

In Figure 6 we plot the phase diagram resulting from the Polyakov lines as a 
function of temperature for different values of the chemical potential, in the fiR-TR 
plane, for = 1, and mR = 0. The line of transitions separates the region in which 
the distribution of the Polyakov line eigenvalues is ungapped from that in which it is 
gapped. The point at which the line of transitions intercepts the Ti?-axis is precisely 
given by the solution of (4.15), which is the result of [14] for = 1, mR = 0. The 
point at which the line of transitions touches the pR-axis is £fi\^j^_Q = 1.5, which is the 
result of [1]. What is perhaps surprising about this result is that the line of transitions 
has very little curvature. This is possible because the line of transitions is not first 
order. Continuity of the effective quark number across the transitions implies that 
they are at least third order. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 6. Phase diagram based on the distribution of the eigenvalues of the Polyakov hne 
for niR = 0, ^ = 1. The line of transitions corresponds to the critical temperatures at 
which a gap forms in the distribution. The quark number is smoothly connected across the 
transitions which implies that they are at least third order. 



It is possible to qualitatively extend the phase diagram of Figure 6 into the region 
of nonzero quark mass. The critical behavior with large mass at zero chemical potential 
is known to approach the Yang-Mills theory result, and at large mass and zero tem- 
perature the critical chemical potential goes to /i = e/i = a/ m'^ + 9/(4/?^), as found in 
[1]. The resulting phase diagram is sketched in Figure 7. 



6 Strong coupling expansion of lattice QCD with /i 7^ 

Making the substitutions (2.6), the Polyakov lines in (5.8) become 

1 ._h{l- x)e-^l^ [8 + 4e2^'/^r2(l + x) - JD(\ - a;)2(3 + x)\ 



1 t\ _ ~ ^y^'^ ^ + 4e-2'^/^r-2(i + x) - JD(1 - xf{?> + x)\ 
~ 16 - JD{1 - x) [16 - JD{1 - x)2(3 + x)] 

where from (5.10) 



4 - JD{1 - x){3 + x) + \/[A- JD{1 - x){3 + x)f - 16(1 - x^h^ 



4(1 -x)e^''^h 



(6.2) 



(6.3) 
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Figure 7. A qualitative picture of the phase transition in fi-T-m space at large N. For /i = 0, 
as m — )• oo the transition goes to the critical temperature for the Yang- Mills theory. For T = 
the lowest transition is set by the energy of lowest fermion mode fi = + 9/(4i?2). There 

are further transitions as [i increases creating a sawtooth pattern and forming lines (shown 

in red) in the T = plane at /u = Jm^^{i^\YlR^, 1=1,2,..., [1]. 



and X is obtained by equating (5.17), in the form 



M ={l + x) 



with (5.21). 



1 Ahr'^e-^"^ [JD{1 - x^) + 2r'^e^^'^{2 - JD{1 - x))] 



1 — X 



+ 



16 - JD{1 - x) [16 - JD{1 - xy{3 + x)] 



(6.4) 



The Polyakov hnes (6.1), (6.2), are plotted in Figure 8 for /i/T = 1.0 and /i/T = 3.0, 
as a function of J. We have set h = 0.01, which corresponds, approximately, to 
ma = 0.05, D = 3, = 4:, and ^ = which we have chosen in the hope of 
comparing to the simulation results in [17-20] (for a review see [21]). However, in the 



very strong coupling limit J = 2 [10] such that the transitions from the strong 

coupling expansion occur for coupling strengths which are too small, {j0-)c — 1-2 for 
yu/T ^ 1.5, while simulation results indicate that the transitions occur for (%r)c — 0.5. 

c 

It would be worthwhile to determine if corrections to the leading contribution in the 
strong coupling and hopping parameter expansions result in a lower {^0-)c- 



Nr 



The effective quark number, = -^Ng, obtained from (6.4) is plotted in Figure 
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Figure 9. Effective quark number M = -^Nq for h = 0.01. (Left) fi/T = 1.0. (Right) 
fj,/T = 3.0. Note that the scale is not the same. 



9 for /i/T = 1.0 and fi/T = 3.0, where the red part of the curve corresponds to the 
ungapped distribution and the blue part to the gapped distribution. It is clear that 
the value of J at which the distribution develops a gap corresponds to the value at 
which the quark number becomes non-zero. Since is always connected during the 
transition from the ungapped to the gapped phase the transitions are at least third 
order. Note that in the lower temperature and larger fi results in [1] it was also found 
that 7^ in the low temperature regions of the ungapped phase when /i > eji, where 
a sawtooth pattern of confinement- deconfinement transitions takes place. 
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0.5 1 1.5 2 2.5 3 3.5 



/i/r 

Figure 10. Phase diagram of large Nc QCD from the lattice strong couphng expansion for 
h = 0.01. The continuity of the effective quark number N implies that the line of transitions 
is at least third order. 

The phase diagram is plotted as a function of J and /i/T in Figure 10, for h = 0.01. 
The deconfinement hne of transitions is defined by the points at which x develops a 
solution which is not —1, or equivalently by the points at which A/" becomes nonzero. 
The line of transitions separates the region in which the distribution of the Polaykov 
hne eigenvalues is ungapped, from that in which it is gapped, and it touches the J axis 
precisely at Jc = -^(1 — 2h) ^ 0.32667, which is the /i = result in [5]. The point at 
which the line of transitions intersects the fi/T axis is given by (yu/T)^ — 3.6. 

Since it is difficult to compare quantitatively with lattice simulation results, even 
when using the same value of h, it is helpful to consider what happens qualitatively as h 
is varied along with J and fi/T, which is sketched in Figure 11. At/i/T = the critical 
behavior is that of [5], J = (1 — 2h), therefore the maximum value of h is |, at which 
the line of transitions as a function of J and n/T reduces to a point. For h = 0, or 
ma = oo, the critical behavior reduces to the Yang- Mills result J = ^. Therefore as 
h is reduced, or as ma is increased, the line of transitions at zero J is stretched along 
the ii/T axis, such that the Yang-Mills result is approached in the /i — )■ limit. 
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h 




/i/T 

Figure 11. A qualitative picture of the phase transitions in fi/T-J-h space. For fi/T = 0, 
the critical line is given by J= -1(1 -2/i) [5], with h 

max — o- For h — the transition is the 
Yang- Mills result J = j^. For J = the transition asymptotes to /i = as /x/T — )• oo. The 
slice of transitions at h = 0.01 as a function of J and fi/T is plotted in Figure 10. 

7 Discussion 

We have calculated the phase diagram of large Nc QCD with small chemical potential, 
^ Sfi, and moderate temperature from one-loop perturbation theory on x 5*^, 
and to leading order in the lattice strong coupling and hopping parameter expansions, 
by solving a single matrix model. It would be interesting to investigate whether it 
is possible to extend both the weak and strong coupling results to the next order, or 
beyond, using a single matrix model and a suitable change of parameters, by including 
also terms with multiply wound Polyakov lines, p±„, for n > 1, which are also known 
to be relevant at higher orders in the strong coupling expansion. Along these lines, 
one could see how the phase diagram is extended into the region of more moderate 
quark masses by including contributions from z^2, Zv3, etc., on x S^, or how the 
phase diagram is extended into the region of larger chemical potentials by including the 
contributions from Zf2, Zfs, etc.. For the lattice theory the former would be equivalent 
to considering additional terms in the strong coupling expansion, and the latter would 
be equivalent to including more terms in the hopping parameter expansion (following, 
for example, [22, 23]). At low temperatures one could expect to reproduce the sawtooth 
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pattern of confinement-deconfinement transitions found in [1]. 

In this work we have considered only the ungapped and single gapped phases where 
the gap opens on the negative real-axis, as expected for the ordinary deconfined phase. 
It might also be interesting to consider multi-gap distributions, which were considered 
in [13] for Yang-Mills theory. 

It might be that approximations exist which lead to a correspondence of matrix 
models at weak and strong coupling for other theories. The qualitative similarity of 
weak-coupling results from 5*^ x 5''^ for QCD with adjoint fermions in [24], and lattice 
simulation results in [25], seems to suggest this possibility. However, in this case it 
could be necessary to consider higher order terms from the strong coupling and/or 
hopping parameter expansion to match the phase diagrams since multi-gap solutions 
are possible, which are a result of transitions in terms with z^n and Zfn for n > 1 in 
the case of the weakly-coupled theory on x S^. 

The behavior of the order of transitions as a function of the chemical potential 
for ordinary QCD with A'^c = 3 is a currently debated issue. If a critical endpoint 
exists at some finite chemical potential where the line of transitions would become first 
order then its location should be bound from below by the radius of convergence of a 
Taylor expansion around fi = of the pressure, the first few terms of which can be 
obtained from lattice simulations, where the radius of convergence is related to the 
quark number susceptibility [26]. However, the predicted location of the lower bound 
on a possible critical endpoint is thought to lie out of the current reach of both lattice 
simulations and experiment [27]. It is interesting that there are simulation results from 
QCD, and the 3-state Potts model, performed in the limit of small chemical potentials, 
which seem to predict that the transition becomes weaker with increasing fj,, suggesting 
that a critical endpoint might not be observed [28, 29]. For these reasons it would be 
interesting to consider further corrections from Zfn ^ for n = 2,3,... to more 
precisely obtain the order of the line of transitions as a function of fi. In Section 3.2.1 
we found that these corrections are only barely perceptible in pi for Yang-Mills theory, 
and that the location and order of the first order phase transition remains unchanged. 
However, it remains possible, and should be checked, if such small corrections to pi are 
able to perturb the Nf ^ system enough to change the order of the phase transition. 
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